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Abstract 



Non - equilibrium tunneling current of an Anderson impurity system subject 
to both constant and alternating electric fields is studied. A time - dependent 
Schrieffer - Wolff transformation maps the time - dependent Anderson Hamil- 
tonian onto a Kondo one. Perturbation expansion in powers of the Kondo 
coupling strength is carried out up to third order, yielding a remarkably sim- 
ple analytical expression for the tunneling current. It is found that the zero - 
bias anomaly is suppressed by the ac - field. Both dc and the first harmonic 
are equally enhanced by the Kondo effect, while the higher harmonics are 
relatively small. These results are shown to be valid also below the Kondo 
temperature. 

The Kondo effect in bulk materials has been extensively studied for more than three 
decades (for review, see [0,0). The impressive advance in the fabrication of mesoscopic 
tunneling devices opens a road to explore the non - equilibrium Kondo physics. Experiments 
have been reported on crossed - wire tungsten junctions ||, quenched lithographic point 
contacts j|, metal and metallic glass break junctions [|J and, recently, quantum dots ||. 
Similar progress has also been recorded in numerous theoretical works, Within this 

realm, theoretical interest is focused on the physics of a Kondo system subject to nonlinear 
time - dependent fields P~12|. Although an appropriate experimental research has not yet 



been carried out, the required ranges of frequencies and temperatures have already been 
achieved in transport measurements so there is no real obstacle to perform the pertinent 
experiments. 

So far, calculations of the current through a Kondo system subject to a time - dependent 
bias were carried out using various assumptions and approximations [0-|TT|]. In the present 
work we carry out straightforward perturbation expansion of the current in powers of the 
coupling constant between the impurity (or the quantum dot) and the conduction bands. 
Indeed, perturbation theory proved its usefulness for the equilibrium Kondo model right 
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from its onset [13]. Recently it was used by Sivan and Wingreen || for the Anderson model 



(in the Kondo regime) with a constant voltage bias. 

Our strategy is to start from the time - dependent Anderson model Hamiltonian and then 
carry out a time - dependent Schrieffer - Wolff transformation to map it onto the correspond- 
ing Kondo - type Hamiltonian. The next step is to combine a specific approach suggested 
by Coleman |14| to treat interaction in strongly correlated systems with the Schwinger - 
Keldysh non - equilibrium Green's functions formalism. It makes the problem amenable for 
perturbation expansion which is performed below up to third order in the Kondo coupling 
J (sixth order in the tunneling strength). The outcome of this procedure is an analytical 
expression for the tunneling current, which is used to obtain two novel results pertaining to 
the direct and alternating currents, namely: 1) The zero - bias anomaly in the direct current 
is smoothed out by an external alternating field. 2) The zeroth and the first harmonic of 
the time - dependent current are enhanced by the Kondo effect while the other harmonics 
remain relatively small. 

At voltages and frequencies less than the level spacing in the tunneling region a convenient 
way to describe a tunneling system is the Anderson Hamiltonian: 

H A = Y ( € k + A L (R){t)) a\ a a k ^ + Y td,aC\ a C d ^ + n ° n v' + Y (ykd<i,a C d,<T + h-C 

keL,R;a cr * rj.a'^u keL,R;a 

(1) 

Here a k a {cLk,a) creates (annihilates) an electron with momentum k and spin o in one of the 
two leads, c da (c dja ) creates (annihilates) an electron with spin a in the tunneling region 
(hereafter we call it "dot"), e k and e dy(T are one - particle energies in the leads and in the 
dot respectively, U is the Coulomb interaction energy in the dot, n a = c da c dtU , while V k d are 
transfer matrix elements between the leads and the dot. The external fields are included 
through the potential shifts of the leads A^^t): 

A L (R) (t) = <j) L{R) + W L{R) cos (ttt + a L ( R )) . (2) 

The first term describes a dc potential, while the second one is due to an ac field, which, for 
simplicity, is assumed to be monochromatic. We note that chemical potentials of the leads 
are shifted by the same amount Al(r) {t) as the one - particle energies, so the population of 
energy levels in the leads does not change. 

Although it is possible to carry out a perturbation expansion within the above model 
using the slave boson technique [§] it turns out to be very cumbersome for time - dependent 
problems. Instead, let us first map the Hamiltonian (11]) onto a Kondo - like one using a 



time - dependent version of the Schrieffer- Wolff transformation ||15|| . This procedure (whose 
details will be explained elsewhere) appears to be straightforward albeit not trivial. Its 
application is limited to the Kondo regime, which is determined by the condition 

e d ,o < 0, e dia + U > 0, \e d>tT \ , e d>a + U >T a , (3) 
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where IV = 2nJ2k^L,R Vkd S(e dtCr — e k ) are the widths of the energy levels in the dot. 
Furthermore, it is assumed that the external fields are not strong enough to draw the system 
out of this regime, so that: 
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4>L(R) 



e d ,a + U\ 



n,W L ,W R < \e d , a \,\e d , a + U\. 



(4) 



It should be stressed however that these conditions do not imply a linear response. The 
later is defined by the conditions 4>l(r) , Wl(r) <C T while T <C \e dj<T \ , \e± a + U\. For the 
sole purpose of having a more compact form of the Hamiltonian the infinite U limit is taken. 
The resulting Hamiltonian then reads, 



H K = 6 k a l,a a k,a + 1/2 Jk'kit) 

k€L,R;a k,k'€L,R;a 



- a a k,<jC dt a C d,-<T + a k' ,a a k,a c d,a C d,a 



(5) 



where 



Jk'kit) = V k > d V k * d exp 
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n 



v £fe' + 0(fc') + s'Q - e djtT e fc + 0( fe ) + sfi - e diCT< 



fed 



fed exp 



sin a(fe) , the symbol (k) means "L" or "R" depending whether 

k belongs to the left or to the right lead, while J B (%) are Bessel's functions. 

Before proceeding with the evaluation of the current, two comments are in order: i) 
Any procedure toward solution of the ensuing transport equations should take into account 
the fact that, out of the full Hilbert space, the system is projected onto a subspace F\ for 
which the dot is occupied by one (and only one) electron, ii) At this stage one might be 
tempted to express the electron creation - annihilation operators in the dot through spin 
operators (T^|, thus arriving at the familiar form |[j,|15| of the Kondo Hamiltonian. But then 
one would realize that the spin operators do not obey the usual commutation rules. In order 
to overcome this obstacle, fictitious (auxiliary) fermions (jTJI might be introduced. But this 
leads one back to equation (0). In other words, auxiliary fermions which sometimes regarded 
as artificial particles introduced to represent spins are real electrons in the dot (impurity 
atom) subject to the constraint specified in i). 

Calculation of tunneling current starting from the Kondo Hamiltonian (|5[) is possible 
for arbitrary field strengths and frequency provided the inequalities @ are satisfied. Yet, 



inspecting a typical experimental setup 
fields and lower frequencies, so that 4>l(r) 



| one may consider somewhat weaker external 
Q, Wl, Wr \ e d ,a\- Expression (]6|) for J k * k then 
significantly simplifies. Matrix elements J k i\~ become time - independent if k and k! belong 
to the same lead. Moreover, J k / k do not depend on potential shifts of every lead separately 
but only on their difference, = Ax — Ar = (p dc + H^cos (Qt + a). For simplicity, we 
further assume that J k i k depend only on the leads to which k! and k belong, independently 
of the values of k' and k. 

The current I(t) is defined as the rate of change in the number of electrons in a lead. 
Within the interaction picture, the commutator of the number operator with the Hamiltonian 
(H) yields an expression for the current, 



J (0 = I Yl Im{j fc /fc(t)Tr Fl [^(4- r J^M(^tW^^(*)5p 



+ 



+ J k , k (t)Tr Fl ^^(a^,^^)^^^)^^)^^)^)]} 



(7) 



(6) 
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where p is the initial (equilibrium) density matrix, while S p and T p are, respectively, the S 
- matrix and the time - ordering operator on a closed time - path. 

We now carry out a perturbation expansion of the above expression in powers of the 
coupling strength using the Schwinger - Keldysh non - equilibrium Green's functions 
technique. In order to get rid of the constraint to the subspace F\ we combine it with the 
method suggested by Coleman [H for the analogous problem in the Anderson model. The 
perturbation diagrams are drawn in Fig. [I]. Every line corresponds to a 2 x 2 matrix of 
Green's functions (details are to be published elsewhere). It is important to note that the 
Hamiltonian ([5]) (and, consequently, the expression (|7|) for the current) contains not only 
the s -d coupling but also the usual potential scattering terms (see Ref. [p]) for explanation). 
As a result, all possible combinations of a, U\ and cr 2 (i.e. o\ = ±o~, cr 2 = ±c) are allowed 
for diagram A but only one of them (namely, a 2 = <J\ = a) is allowed for diagram B. The 
spin - flip diagrams of the type B cancel. 

The resulting current is given by the following formulae: 



I(t) =I<V(t) +I {3) (t) 
j( 2 )(t) = C 2 [<p dc + W cos(fit + a) 

7 (3) (t) = -To + Yj l J «l cos(nfii + na + arg I n ), 

Z n=l 



(8) 



+oo 



^ ,w w 

C3 2^ Js(-^)Js+n(-^) 



idc 



+ sfi, T, D) 



n 77* / idc 



-ITF 



sn,T,D] 



F(<P,T,D) = -Re 



duo de 



[f L (co)-f R (cu)}[f L (e) + f R (e)] . 



D 



uo — e + 27 



in — coth — 
2 2 



where 1^ and express contributions of second (diagram C) and third (diagrams A and B) 
orders in J yk p, while C 2 = |7T J lr p L p R and C 3 = §57r J lr p L p R (JllPl + JrrPr)- The 



quantities pi{R) are densities of states in the leads, whereas Jl(r)r(l) = VL(R)V R (L)/\ e d,a\- The 
cutoff D is equal to the energy difference between the chemical potential and the bottom of 
the conduction band, while /l(c) = 1/ (exp[(e — 4>)/kT] + 1) and f R (e) = 1/ (exp[e/A;T] + 1) 
are Fermi functions in the leads (the left lead being shifted by 0). The integral in the 
expression for F can be written in the same form as in the dc result of Sivan and Wingreen || . 
The alternating field causes splitting of the energy levels in one of the leads [0] . Therefore 
the time - dependent current is a result of interference between "dc - like" contributions, 
each one with an effective bias (f) dc ± sfl. 

An approximate evaluation of the double integral in the above formulae is possible both 
for the linear (<p T) and for the non - linear (0 3> T) regimes. The result is 



F(<P,T,D) 



if 



,ln f' 



< T 
> T 



(9) 



The terms that are not included in this expression are of order <p. It should be stressed that 
in the nonlinear case on which we focus our attention here the function F does not diverge 
as T —>■ 0. Hence, our results for the nonlinear response must be valid even below the Kondo 
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temperature. This results from the fact that here, the non - linear bias plays the role of 
temperature as the largest low - energy scale. 

Substitution of expression (^) into equations (Q) yields simple formulae for the current. 
First of all, in the absence of an ac - field we obtain, 



C 2 + C 3 \n§, if dc «T, 
C 2 + C 3 \nj^, if dc >T 



(10) 

I0 dc " 



We note that, to the best of our knowledge, such a simple expression for a non - equilibrium 
tunneling current through a Kondo system has not been found before. In the presence of a 
strong (W ^> fl) ac - field we get, 

I{t) « I dc + I ac {t) (11) 

C 2 + C 3 In- ° 



rdc idc 



A dc (^ } W,n } T) 

I ac (t) = Wcos(nt + a) 



C 2 + C 3 \n D 



A ac (<^,w,n,T) 

where Aa c , A ac are of the order of max(0 dc , W), if <ft dc > T or W > T (nonlinear response), 
and A dc = A ac = T, if (p dc <C T and IfcT (linear response). 

In most of the relevant experiments the attention was focused on the dc. The best 
pronounced feature of the Kondo effect in the dc seems to be the zero - bias anomaly, i.e. 
appearance of a narrow peak in the differential conductance around zero bias. In Fig. |2] and 
[3] we show how it is altered by application of the external time - dependent field. According 
to formula (§), dc in an external ac - field is given by a sum of direct currents without ac - 
field but with effective biases (p dc + sQ (these currents are weighted by Jf). Consequently, 
the zero - bias anomaly is suppressed. A picture of gradual flattening similar to Fig. |2| holds 
also at lower frequencies down to the adiabatic limit. At higher frequencies, however, side 
peaks at multiples of Q can be resolved (Fig. |3|). 

Although measurement of an alternating current with frequencies and amplitudes in the 
relevant range is not an easy task, it might reveal new interesting features of the Kondo 
effect. Expression (|TT|) shows that only the dc and the first harmonic are enhanced by 
the Kondo effect. For the other harmonics, the interference of the contributions to the 
current in equations @ with different effective biases is destructive. We emphasize that 
the nonlinear response is different from the linear one. Namely, (1) the second and the 
higher harmonics exist although they are not amplified by the large factor ln-^, (2) the 
factor In ^ is replaced by ln-^. Nevertheless, they look very similar in the sense that the 
dc and the first harmonic are much larger than the others. This is a particular feature of 
the interacting (Kondo) system. In the non - interacting one - level system the amplitude 
of the higher harmonics I n , n > 2 is, generally, of the order of the amplitudes of the dc 
and the first one (gradually decreasing with their number as V / (nfl) ) . In some sense the 
Kondo system, in contrast with the non - interacting one, behaves like an ordinary resistor 
(although the current is enhanced by the Kondo effect). Unlike the situation encountered 
in non - interacting systems, the magnitudes of the dc and the ac components are mostly 
governed by two independent parameters cj) dc and W . Experimentally, Kondo contribution 
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to the tunneling current is usually revealed through a special dependence on the parameters 
(such as InT increase of the conductance or zero - bias anomaly). Our formulae (]TTf) imply 
that this kind of effects can be found in the first harmonic as well as in the dc but not in 
higher harmonics. 

In conclusion, the pertinent physics is rather rich since it combines strong correlations 
with nonlinearity and time dependence. We obtain (within perturbation theory) rigorous an- 
alytical formulae for the nonlinear time - dependent tunneling current in the Kondo regime. 
They apply for nonlinear response both below and above the Kondo temperature (for linear 
response they are valid only above it). In the special case of constant external field the I — (f) 
relation ( PD encodes the main physics of the familiar zero - bias anomaly. 

The first novel result of the present research can be easily tested experimentally, since 
it concerns the direct tunneling current. We find that its zero - bias anomaly is suppressed 
by an ac - field. Second, it is shown that both the zeroth and the first harmonics of the 
alternating current are strongly enhanced by the Kondo effect, while the other harmonics 
remain relatively small. This result is slightly more difficult to be tested experimentally. 
Yet, it is apparently worth the effort since it is remarkably different from what is expected 
for a non - interacting one - level system where all the harmonics emerge together. 

As far as relation to previous relevant works is concerned, we, first, notice that our ana- 
lytical results for the dc are consistent with the numerical calculations of Ref . || . However, 
being able to consider stronger ac - fields (larger ratio W/Q), we find also an overall suppres- 
sion of the zero - bias anomaly, besides the appearance of side peaks. As for the spectrum 
of the tunneling current, we are unable to validate the assumption suggested in Ref. || that 
all the harmonics beside the dc one can be neglected. On the other hand, we verify that the 
second and higher harmonics are generated but they are indeed much smaller than the dc 
and the first one [P|. Within a specific model, some authors |TU[ obtained current spectrum 



similar to that of a non - interacting system. We attribute the difference between this result 
and ours to the very peculiar choice of parameters used therein. 
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FIGURES 



FIG. 1. Diagrams for the perturbation expansion of the current (^). Solid lines stand for lead 

electrons, dashed lines — for dot electrons (auxiliary fermions). 

FIG. 2. Suppression of the zero - bias anomaly in the direct current by an external alternating 

field at relatively low frequencies. Differential conductance dl dc /d(f) dc is given in units of C3T, 

while D = 50, HQ, = 5, W and (fi dc are in units of T. 

FIG. 3. Suppression of the zero - bias anomaly in the direct current by an external alternating 

field at higher frequencies. Side peaks appear at multiples of f2. Differential conductance dl dc / d<j) dc 

is given in units of C3T, while D = 200, HQ = 15, W and <j) dc are in units of T. 



